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ABSTRACT 

Using pulse timing observations we have obtained precise parameters, including positions with ~ 20 mas ac- 
curacy, of five millisecond pulsars in NGC 6752. Three of them, located relatively close to the cluster center, 
have line-of-sight accelerations larger than the maximum value predicted by the central mass density derived from 
optical observation, providing dynamical evidence for a central mass-to-light ratio > 10, much higher than for 
any other globular cluster. It is likely that the other two millisecond pulsars have been ejected out of the core to 
their present locations at 1.4 and 3.3 half-mass radii, respectively, suggesting unusual non-thermal dynamics in 
the cluster core. 

Subject headings: Globular clusters: individual (NGC 6752) — pulsars: individual (J1911-5958A; J1910-5959B; 
J191 1-6000C; J1910-5959D; J1910-5959E) 



1. INTRODUCTION 

During a deep search of the globular cluster (GC) system for 
millisecond pulsars (MSPs), carried out using the Parkes ra- 
dio telescope, we discovered (D'Amico et al. 2001a) a binary 
MSP (hereafter PSR-A) associated with NGC 6752, a core- 
collapsed cluster showing evidence of mass segregation (Fer- 
raro et al. 1997). The resulting knowledge of the dispersion 
measure (DM) for this cluster has facilitated the discovery of 
four additional MSPs (hereafter PSRs B, C, D, E) in the same 
cluster (D'Amico et al. 2001b). We have made frequent ob- 
servations using the Parkes telescope to obtain coherent timing 
solutions for these pulsars. These results can be used to esti- 
mate a variety of physical properties of the pulsars and of the 
host cluster. In particular, if the measured period derivatives are 
dominated by the dynamical effects of the cluster gravitational 
field, they can be used to constrain the mass-density distribution 
of the cluster, giving information on the GC's dynamical status 
and on the population of optically unseen cluster members. In 
this paper we report on the first 15 months of pulse timing ob- 
servations, providing rotational and positional parameters for 
the five MSPs known in NGC 6752, and we use the pulsar po- 
sitions and inferred accelerations to probe the dynamical status 
of the cluster. 

2. OBSERVATIONS AND RESULTS 

Regular pulsar timing observations of NGC 6752 have been 
carried out since 2000 September with the Parkes 64-m radio 
telescope, using the center beam of the multibeam receiver at 
1400 MHz. The hardware system is the same as that used in the 
discovery observations (D'Amico et al. 2001a). The effects of 
interstellar dispersion are minimized by using a filter bank hav- 
ing 512 x 0.5 MHz frequency channels for each polarization. 
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After detection, the signals from individual channels are added 
in polarization pairs, integrated, 1 bit-digitized every 125 /is, 
and recorded to magnetic tape for off-line analysis. Observation 
times are typically 1-2 hours. Pulse times of arrival (TOAs) are 
determined by fitting a standard high signal-to-noise ratio pulse 
profile to the observed mean pulse profiles and analyzed using 
the program TEMPO 1 and the DE200 solar system ephemeris 
(Standish 1982). 

Table 1 lists the timing parameters obtained for the five 
pulsars, including precise positions. Because these pulsars 
have relatively low DMs^33 cm" 3 pc, interstellar scintillation 
strongly affects their detectability. In particular, PSR-B was 
detected on only 27 of 140 observations of the cluster. The 
mean flux densities at 1400 MHz (S1400) in Table 1 are aver- 
aged values, derived from the observed signal-to-noise ratios. 
Non-detections were accounted for by assuming S1400 values 
corresponding to half the detection limits. As can be seen from 
Table 1, all but PSR-A are isolated MSPs. PSRs B, D and E are 
located close to the cluster center. PSR-D has the third largest 
period derivative, P = 9.6 x 10~ 19 , among known MSPs, after 
PSR B1820-30A in NGC 6624 and PSR B 182 1-24 in M28, 
suggesting that the P is dominated by the line-of-sight accel- 
eration in the cluster gravitational field. This interpretation is 
supported by the large negative P values observed for PSR B 
and E, which are also located close to the cluster center. 

PSR-C is located 2'.1 from the cluster center, equivalent to 
about 1 .4 half-mass radii or 24 core-radii, assuming a core ra- 
dius 9 C = 6."7 (Lugger, Cohn & Grindlay 1995), a half-mass 
radius 9i„„ = 115" (Djorgovski 1993), and an optical center for 
the cluster of (J2000) R.A.=19 h 10 m 51?8 and Dec.=-59° 58' 
55" (Harris 1996). Previously, the largest offset of an associ- 
ated pulsar from a GC center was for PSR B2 127+1 1C (Prince 
et al. 1 99 1 ), a member of a double-neutron-star eccentric binary 
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in M15, located at ~ 0.8 half-mass radii. Given the large offset 
of PSR-C from the GC center, its period derivative is not sig- 
nificantly affected by the GC potential well, allowing the mea- 
surement of the characteristic age (r c = P/2P = 3.8 x 10 10 yr, 
the largest among known pulsars), the surface magnetic field 
(B sm£ = 3.2 x 10 19 (P^) 1/2 = 1.1 x 10 s Gauss) and the rotational 
energy loss (E = 3.95 x 10 46 P/f 3 =5.9 x 10 32 erg s" 1 ). 

PSR-A, the first pulsar discovered in this cluster, is located 
even farther from the center, at 6'4, equivalent to 3.3 half-mass 
radii or ~ 57 core radii. Given this large radial offset, one could 
question the association of the pulsar with the cluster. Based on 
the 19 MSPs with S 436 £ 2 mJy detected by the Parkes Southern 
Pulsar Survey (Lyne et al. 1998) and assuming a typical spec- 
tral index ~ -1.9 (Toscano et al. 1998) for MSPs, the prob- 
ability of chance superposition of a Galactic field MSP with 
5i4oo ~ 0.2 mJy within 6.'4 of the center of NGC 6752 is at 
most 10~ 5 . This probability is further reduced by noting that 
all five pulsars have very similar DM values. Also in this case 
the observed period derivative can be used to measure the pul- 
sar parameters r c = 1.7 x 10 10 yr, B sur f = 1.0 x 10 8 Gauss, and 
£ = 5.9 x 10 32 ergs" 1 . 

Accurate DMs were obtained for each pulsar by splitting the 
256-MHz bandwidth into two adjacent sub-bands and comput- 
ing the differential delays. The maximum deviation (~0.3 cm" 3 
pc) of the DM values derived for each pulsar from the aver- 
age (33.36 cm" 3 pc) is similar to that observed in other GCs 
hosting several MSPs like 47 Tucanae (Freire et al. 2001) and 
M15 (Anderson et al. 1990). Given the wide angular offsets 
of some of the pulsars, the scatter in some of the DMs could 
arise from gradients in the Galactic column density across dif- 
ferent lines-of-sight toward the cluster (Armstrong, Rickett & 
Spangler 1995). It could also arise from an enhanced electron 
density within the cluster (see Freire et al. 2001). Excluding 
the peripheral PSR-A, a rough estimate of the density of such 
gas is given by n e ~ ((ADM) 2 } 1 / 2 /(D(6» 2 L ) 1 / 2 ) = 0.025 ±0.005 
cm , where ((ADM) 2 ) 1 / 2 is the rms deviation in the DMs and 
{6\) l l 2 is the one-dimensional dispersion of the angular offsets 
in radians with respect to the GC center in the plane of the sky 
for the 4 inner MSPs. D = 4.1 ±0.1 kpc (Renzini et al. 1996) 
is the cluster distance (the errors are reported at the lc-level, as 
everywhere in this paper and in Table 1). This n e value is about 
a factor of ten smaller than that inferred from the four MSPs 
(with P < 10 ms) in the core-collapsed cluster M15 and ~ 40% 
of that derived using a more refined modeling of the plasma 
content in 47 Tuc (Freire et al. 2001). 

NGC 6752 was observed recently for ~ 28700 s with the 
ACIS-S detector aboard the Chandra X-ray observatory by 
Pooley et al. (2002). The position of PSR-D is consistent with 
that of the Chandra source labeled CX11 by Pooley et al. All 
the MSPs but PSR-A are located in the Chandra field of view, 
but PSR-C is outside the half-mass radius region searched for 
X-ray sources by Pooley et al. Thus we have processed the full 
ACIS-S3 image in the 0.5-6.0 keV band using the CIAO 2.2 
software 2 , resulting in the detection of two additional proba- 
ble X-ray counterparts to the MSPs. Using the wavedetect tool, 
we have found a weak soft source whose error circle encloses 
PSR-C, having a hardness ratio (as defined by Grindlay et al. 
2001a) > 5.5. Assuming isotropic X-ray emission and a hydro- 
gen column density Nh = 2.2 x 10 20 cm" 2 , we obtain 3 for three 
different spectral models (a power-law of photon index -2.5, 



a blackbody with kT = 0.3 keV or a thermal bremsstrahlung 
with kT = 1 keV) similar values of the X-ray luminosity in the 
0.5-2.5 keV band, Lx=2.2xl0 30 erg s" 1 , corresponding to a 
conversion efficiency Lx/E ~ 0.004. This is somewhat higher 
than that predicted on the basis of both the sample of MSPs 
observed in 47 Tuc and the sample of MSPs in the Galactic 
field (Grindlay et al. 2002). Using the celldetect tool, we have 
also found marginal evidence for a slightly harder source com- 
patible with the position of PSR-B, having a hardness ratio ~ 
4 and Lx = 1.1 xlO 30 erg s" 1 . This source is surrounded by 
many other sources and the nominal celldetect signal-to-noise 
ratio, ~1.5, is low. No X-ray source is associated with PSR- 
E, with an upper limit to the X-ray luminosity in the 0.5-2.5 
keV band of 10 30 erg s . Based on the three probable posi- 
tional coincidences, we have used the MSP positions obtained 
via radio timing to derive a corrected Chandra astrometric solu- 
tion, requiring a shift of the Chandra positions by -0"05 in RA 
and -0."03 in Dec, with a final rms (of the differences between 
Chandra and radio positions) of 0."28 in RA and 0"39 in Dec, 
consistent with the expected positional uncertainties of ~ 0"3 
(Pooley et al. 2002). 

Figure 1 shows an optical image of the cluster, the positions 
of the five MSPs and the contour levels of the Chandra im- 
age. Pooley et al. (2002) pointed out that the X-ray sources 
in the central 1 .'6 x 1 .'6 region lie in the south-east quadrant. 
While they ascribe this to chance, we note that most of the X- 
ray sources within the half-mass radius region are roughly dis- 
tributed along a stretched S-shaped pattern, whose elongated 
ends are oriented in the east-west direction. 



3. CENTRAL MASS-TO-LIGHT RATIO 

The large negative P values observed in PSRs B and E can 
be used to derive lower limits to the line-of-sight accelerations, 
fl ,/c = P/P = -9.6±0.1 x 10" 17 s" 1 , for both pulsars, which are 
the largest known after those of PSRs B2 127+11 A and D in 
the core of M15 (Anderson et al. 1990). Contributions from 
centrifugal acceleration (Shklovskii 1970), differential Galac- 
tic rotation (Damour & Taylor 1991) and vertical acceleration 
in the Galactic potential (Kuijken & Gilmore 1989) are all neg- 
ligible. Finally, according to Figures 3 and 4 of Phinney ( 1 993), 
the probability that the accelerations of both PSRs B and E are 
dominated by a nearby star in the cluster is < 10~ 4 . Thus we 
can conclude that the inferred high values of \ai/c\ are due to 
the potential well of NGC 6752. 

A lower limit to the mass-to-light ratio in the inner regions of 
NGC 6752 can be derived from the following rule, which holds 
to within ^10% in all plausible cluster models (Phinney 1992): 
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Here Sy(< 8±) is the mean surface brightness within a line of 
sight subtended by an angle 6± with respect to the cluster cen- 
ter, M cy i(< 8±) is the mass enclosed in the cylindrical volume 
of radius R± = D8± and M. / £\> is the mean projected mass-to- 
light ratio in the V-band. 
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Fig. 1. — Positions of the 5 millisecond pulsars (arrows) and of the 
X-ray sources (black contours) in NGC 6752, superimposed on an op- 
tical image retrieved from the Digital Sky Survey. The circle indicates 
the half-mass radius region (Ohm = 115"). The rectangular box indi- 
cates the Chandra ACIS-S3 field of view. North is to the top and east 
is to the left. 

We use the most recent published brightness profile for this 
cluster (Lugger, Cohn & Grindlay 1995), normalized to the cen- 
tral surface brightness in the V-band reported by Djorgovski 
(1993), in order to plot the curves of maximum \a//c\ for differ- 
ent values of M.j C\>. In Figure 2 the histogram represents the 
data for M/ '£y=l.l, as suggested by Pryor & Meylan (1993). 
This greatly underestimates the observed values of P/P. The 
dashed lines are analytical fits (with reduced \ 2 ~1) to the ob- 
served data, scaled according to increasing values of M./ C\>. 
Only M /C v >9 can account for the observed \P/P\ of PSRs B 
and E. An even larger M./ C\> ^ 13 is required if PSR-D has a 
negligible intrinsic P/P, as might be expected given the close 
projected positions and the similar absolute values of P/P for 
PSRs B, D and E. A reasonable interpretation is that PSR D is 
at about the same distance from the GC center as the PSRs B 
and E, but in the closer half of the cluster, whereas PSRs B and 
E reside in the further half. In this case, for these 3 MSPs the in- 
trinsic P/P ~ 5 x 10~ 18 can be estimated as the average of those 
of PSRs B (or E) and D. Values of M / C v in the interval 9-13 
result if the observed scalings between X-ray luminosity and 
spin-down power for MSPs are used to estimate the magnetic- 
dipole braking contribution to P/P for PSR-D (see caption of 
Fig. 2). We finally note that these estimates of the projected 
M. / £y are independent of distance, excepting the effects of 
extinction, which are very small for NGC 6752, with E(B-V) = 
0.04 (Harris 1996), and of modeling of the cluster potential. 

These results imply (see eq. |l]) that there must be M cy i(< 
0±,e) ~ 1-3 x 10 4 M Q of matter in the form of low-luminosity 
stellar objects within the projected radius of PSR-E, 0.15 pc. 
Adopting the prescription of Djorgovski (1993) and a core ra- 
dius 9 C = b'.'l (Lugger, Cohn & Grindlay 1995) this in turn cor- 
responds to a central mass density J7x 10 5 M Q pc" 3 , at least 
five times larger than that derived from measurements in the op- 
tical band (Pryor & Meylan 1993). The stars whose initial mass 
was in the interval 0.6-0.8 M Q (bracketing the current turn-off 
point of the cluster) now dominate the total integrated V-band 
luminosity of NGC 6752, 1.2 x 10 5 L V0 (Djorgovski 1993). 



M/L v = 3 




M/L v =1.1- 



5 10 15 20 

8 [ arcsec ] 

Fig. 2. — Maximum line-of-sight acceleration \ai ma Jc\ = \P/P\ ver- 
sus displacement with respect to the center of NGC 6752. The his- 
togram represents the prediction based on the available optical obser- 
vations. The dashed lines are analytical fits to the optical observations, 
labeled according to the adopted mass-to-light ratio. The filled pen- 
tagons represent lower limits to the line-of-sight acceleration based on 
the observation of PSRs B and E. The open pentagon shows the value 
of P/P for PSR-D assuming a negligible intrinsic P. If various ob- 
served scalings (Becker & Triimper 1997; Possenti et al. 2002; Grind- 
lay et al. 2002) between X-ray luminosity and spin-down power for 
the MSPs are used to estimate the intrinsic P, P/P for PSR-D assumes 
values in the interval between the open and the dotted pentagon (see 
§3). 

Assuming a typical mass-to-light ratio for these stars of 
~ 0.25 (Phinney 1993), the number of 0.6-0.8 M stars is 
N t0 ~4x 10 4 . Assuming a Salpeter initial mass function (IMF) 
oc m~ a with a = 2.35, the evolution of stars having initial mass 
Z 8M Q should have left at most Ndr ~ 4 x 10 3 dark remnants. If 
they sank towards the cluster center and are responsible for the 
high value of the mass-to-light ratio in the core, their average 
mass would be (m^) =M cy i(< 6±.E)/Ndr Z 3.2 M Q , suggesting 
that at least some of them must be black holes. The same con- 
clusions apply to any power-law IMF steeper than Salpeter's. 
Conversely, if we assume that the low-luminosity objects within 
the central 0.15 pc are all 1.4M Q neutron stars (NSs), this 
would imply «ns < 2.1 for the slope of the IMF. Given the 
very low upper limit on the density of any ionized gas, > 10 4 
NSs could reside in the inner region of NGC 6752 without pro- 
ducing an accretion luminosity detectable in the deep Chandra 
exposure of Pooley et al. (2002). Alternatively, if the popula- 
tion of the central region is dominated by massive white dwarfs 
(descendents of stars of initial mass Z 2.5 M Q ), we obtain a 
nominal limit «wd < 2.5. However, these ~1M Q remnants 
(less massive than the population of MSPs) should distribute 
out to larger distances from the cluster center than PSR-E. A 
conservative 60% increase in their estimated number implies 
an IMF flatter than Salpeter's in this case as well. 

We note that both the uncertainties in the average mass-to- 
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light ratio of the 0.6-0.8 M Q stars and the unknown number of 
these stars which have been tidally stripped or evaporated from 
the cluster since its formation could affect N to and thus our re- 
sults. However, simulations of strongly concentrated clusters 
(e.g. Joshi, Nave & Rasio 2001) indicate that the relative deple- 
tion of 0.6-0.8 M Q stars should be smaller than that of neutron 
stars, which escape the cluster potential well due to supernova 
kicks at birth (Rappaport et al. 2001). Such an effect would 
further increase the average mass (mar) of dark remnant stars 
and/or make the IMF even flatter. 



4. DISCUSSION 

Our results provide the first direct dynamical evidence for a 
high density of unseen remnants in the core of a globular clus- 
ter, suggesting a mass-to-light ratio > 10 in the core region. For 
comparison, we note that a much smaller value ~ 3, was ob- 
tained for the case of M15 (Phinney 1993), a core-collapsed 
GC long suspected to host a central black-hole (see e.g. van 
der Marel & Roeland 1999). The nature of the remnants in the 
core of NGC 6752 is not clear as it strongly depends on as- 
sumptions about the IMF. However, the possibility that they are 
black holes, or that many stellar remnants have collapsed into a 
single massive black hole, is intriguing. 

In addition, the large offsets of PSRs A and C from the cluster 
center indicate the occurrence of highly effective non-thermal 
dynamics in the cluster core. No other GC shows an MSP 
ejected beyond its half-mass radius (which requires a finely 
tuned impulse to avoid prompt expulsion of the NS from the 
GC), while NGC 6752 has two of them. They could result from 
exchange encounters in the core (Phinney & Sigurdsson 1991) 
but the large offset of PSR-A suggests also the occurrence of 



more powerful scattering events (Colpi, Possenti & Gualandris 
2002). The scattering target must have treated this binary sys- 
tem as a point mass, ejecting it without inducing appreciable 
eccentricity. Assuming a value of 1.4 M Q for the pulsar, the 
total mass of the binary system containing PSR-A is at least 1 .6 
M Q . Simple dynamical considerations favor a scattering target 
significantly more massive than the scattered binary, supporting 
the conclusion that it could be of many solar masses. A black- 
hole binary system would be a natural candidate. The Chandra 
X-ray observations do not place severe limits on its maximum 
mass. The arguments of Grindlay et al. (2001a) applied to an 
electron gas density < 0.025 cm" 3 and a detection threshold of 
~ 10 30 erg s" 1 , allow for up to a few hundred solar masses in the 
form of black hole(s) present in the central region of NGC 6752. 

Finally, the pattern of X-ray sources, roughly oriented mid- 
way between the probable projected ejection directions of PSRs 
A and C, suggests a preferential geometry for ejection, similar 
to that discussed by Grindlay et al. (2001b) for NGC 6397. Its 
physical connection with the inferred high density of unseen 
matter in the core is not clear, but it could be a further signature 
of significant non-thermal activity in the core region. 
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ica Colpi, and Piero Ranalli, Marcella Brusa and Paolo Montegriffo for 
assistance with reduction of the X-ray data and astrometry. N.D' A. and 
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and the Italian Minister of Research (MIUR). F.C. acknowledges sup- 
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Table 1 

Measured and Derived Parameters for Five Millisecond Pulsars in NGC 6752 



Pulsar 


PSR-A 


PSR-B 


PSR-C 


PSR-D 


PSR-E 


Name 


J1911-5958A 


J1910-5959B 


J1911-6000C 


J1910-5959D 


J1910-5959E 


R.A. (J2000) 


19 h ll m 42. s 7562(2) 


19 h 10 m 52?050(4) 19 h ll m 05?5561(7) 19 h 10 m 52?417(2) 19 h 10 m 52?155(2) 


Dec (J2000) 


-59° 58'26."900(2) 


-59° 59'00."83(3) 


-60° 00'59"680(7) 


-59° 59'05"45(2) 


-59° 59'02."09(2) 


P(ms) 


3.2661865707911(5) 


8.35779850080(3) 5.277326932317(4) 


9.03528524779(2) 4.571765939765(7) 


P 


3.07(10)xl0" 21 


-7.99(5) xlO" 19 


2.2(7) xlO -21 


9.63(3) xlO" 19 


-4.37(l)xl0" 19 


Epoch (MJD) 


51920.0 


52000.0 


51910.0 


51910.0 


51910.0 


DM (cm~ 3 pc) 


33.68(1) 


33.28(4) 


33.21(4) 


33.32(5) 


33.29(5) 


Porb (d) 


0.837113476(1) 










a v sin(/) jc (s) 


1.206045(2) 










T asc (MJD) 


51919.2064780(3) 










Eccentricity 


< 10" 5 










M c (M Q ) 


> 0.19 










MJD Range 


51710-52200 


51745-52202 


51710-52201 


51744-52197 


51744-52201 


No. TOAs 


74 


27 


94 


38 


38 


residual (/is) 


10 


83 


55 


55 


60 


5i4oo (mJy) 
Offset a (') 


0.22 


0.06 


0.30 


0.07 


0.09 


6. '39 


O.'IO 


2f70 


0.'19 
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Angular separation in the plane of the sky between the MSP and the center of NGC 6752 (Harris 1996). 



